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NEW IDENTITIES INVOLVING 
BERNOULLI AND EULER POLYNOMIALS 

Hao Pan and Zhi-Wei Sun 



Abstract. In this paper we obtain several new identities for Bernoulli 
and Euler polynomials; some of them extend Miki's and Matiyasevich's 
I identities. Our new method involves differences and derivatives of polyno- 

mials. 



^ ■ 1. Introduction 



The Bernoulli numbers Bq, Bi, . . . are given by 



(N 

> 

cn 



1 = ^^^^;^ (N<2^) 



they play important roles in many aspects. Most research on Bernoulli 
O . numbers concentrates on their congruence properties (cf. e.g. [Sul]). 
. However, there are also some interesting identities concerning Bernoulli 
numbers (see, e.g. [Di] and [Su2]). 
15 . In 1978 Miki [Mi] proposed the following curious identity which involves 

both an ordinary convolution and a binomial convolution of Bernoulli num- 
bers: 



BhB„-u \ ^ / n\ BiBn-i 2 



for any n = 4, 5, . . . , where 



1 1 

i^n = l + - + ■■•+-. 

2 n 

In the original proof of this identity, Miki showed that the two sides of 
(1.1) are congruent modulo all sufficiently large primes. Shiratani and 
Yokoyama [SY] gave another proof of (1.1) by p-adic analysis, and Gessel 
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[Ge] reproved Miki's identity (1.1) by using the ordinary generating func- 
tion and the exponential generating function of Stirhng numbers of the 
second kind. 

Inspired by Miki's work, Matiyasevich [Ma] found the following two 
identities of the same nature by the software Mathematica. 

—j-Bn-k - ^ f ^ j -yBn-l = HnBn (1.2) 
k=2 1=2 ^ ^ 

and 

n-2 n-2 y -L 2\ 

(n + 2) ^ BkBn-k -^J2r 1 ) = + 1)^" (1-^) 

k=2 1=2 ^ ^ 

for each n = 4,5, We mention that (1.2) is actually equivalent to 

Miki's identity (1.1). The reason is as follows: 

BkBn-k _ \ ^ / n\ BiBn-i 

k=2 ^ ^ 1=1 \ \ 

2 ^ Sfe ^ fn\Bi 



= nT.l^Bn-k--J2[l)T^n-l. 
k=2 1=2 ^ ^ 

Quite recently Dunne and Schubert [DS] presented a new approach to 
(1.1) and (1.3) motivated by quantum field theory and string theory. 

The BernouUi polynomials Bn{x) = ^Lo (fc)^^^""'' (n = 0, 1, 2, . . . ) 
have the following basic properties: 

Br,{x + 1) - Bn{x) = nx^-\ Bn{x + y) = J2 Cf)Bk{x)y''-'', 

k=0 ^ ^ 

and also B'^{x) = nBn-i{x) for n > 0. 

In this paper we extend Miki's identity (1.1) and Matiyasevich's identity 
(1.3) to Bernoulli polynomials. 

Theorem 1.1. Let n > 1 be an integer. Then 

Bk{x)Bn-k{y) - ^\ Bijx - y)Bn-i{y) + Bi{y - x)Bn-i{x) 

= ^(i?.(x)+i?.(,)) + ^-(f-^f) 
n n{x — y) 

(1.4) 
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and 

'n + 1\ Bi{x - y)Bn-i{y) + Bi{y - x)B 



Y.Bk{x)Br.-M-Y.\^^ll) 



1 + 2 

^ Bn+l{x) + Bn+l{y) _ 2 _ Bn+2{x) - -gn+2(j/) 

(x — y)'^ n + 2 (x — y)^ 

(1.5) 

We remark that (1.4) has the foUowing equivalent version 



n-l 



Bk{x)^ ^ ^ N , D ^ Bn{x)-B^{y) 



2 E ^^Bn-k{y) - Hn-i{B^{x) + BM) 



Corollary 1.1. Le^ 2 he an integer. Then we have 



fe=l ^ «=2 ^ ^ 

and 



Y^B,{x)B^_,{x)-2j2 (?+i)^TTy^ = (n+l)S„(a;). (1.7) 

Proof. Letting y tend to x and recaUing that B'^{x) = nBn-i{x), we 
immediately find that (1.4) implies (1.6). 

Now we go to prove (1.7). Let P{z) = Bn+2{z)/{n + 2). Then P'{z) = 
Br,+i{z), P"{z) = {n+l)Bn{z) and P"'{z) = n{n + l)Bn-i{z). In light 
of Taylor's expansion, 

P{y) - P{x) = P'{x){y -x) + -^{y - xf + -^{y - xf + ■ ■ ■ 



and 



P\y) - P'{x) = P"{x){y - x) + -^{y - xf + 
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Therefore 

■ Bn+l {x) + Bn+l {y) 2 Bn+2 (x) - Bn+2 (v) 



lim 



— lim 



y- 



{x — n + 2 {x — yY 

P'{x) + P'{y) 2{P{x) ~ P{y)) \ 



x\ {x-yY {x-yY J 



y^x \ (y — xY {y — xY 

fP"{x) P"'ix) ^( P"{x) P"'ix) 

= lim ^ + + 2 ^ \ + + 

y^x\y-x 2! \2l{y-x) 3! 



6 6 

In view of this, we can easily get (1.7) from (1.5) by letting y tend to 

X. □ 

Now let us see how Miki's identity follows from (1.6). In fact, (1.6) in 
the case a; = yields that 

EBuB„—k \ — ^ f — 1\ BiBr,—] 2B„ „ „ 2 „ 

j^-i /c(n — k) ^ \ I J l{n — l)n'^ n 

1=1 

n-l 



/ J \n — I J J l{n — I) n 

- £ (;) 1^ + l^nBn + 5.-1. 



\Zy l{n — I) n 

We also mention that (1.7) in the — gives Matiyasevich's identity 

(1.3). 

Corollary 1.2. Let n ^ A be an integer. Then 

E"" ^ Bk ^ _ n ^ BkBn-k ^ fn\Bk ^ ^ 

fe=2 fc=2 ^ fc=2 / 

where Bk = (2^"^ - l)Bk. 

Proof. Simply take x = 1/2 m (1.6) and use the known formula S.(l/2) = 
B^. (Note also that n/{k{n - k)) = 1/k + l/{n - k).) □ 

The last equality in Corollary 1.2 was first found by C. Faber and R. 
Pandharipande, and then confirmed by Zagier (cf. [FP]). 
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The Euler polynomials En{x) (n = 0, 1, 2, . . . ) are defined by 



n=0 



Here are some basic properties of Euler polynomials: 

and also E'^{x) — nEn-i{x) if n > 0. It is also known that 

En{x) ^^(^B^+i{x)-2-+'B^+, (1)^. 

Similar to Theorem 1.1 we have the following identities involving Euler 
polynomials. 

Theorem 1.2. Let n be a positive integer. Then 



n + 2 X — y 



n+l 



„^(n+ l\ E,(x - y)B„+i_,(y) + 1 



_i{y) + Ei{y- x)Bn+i-i{x) 



(1.8) 



Also, 



1=1 

and 



± ^E,.M - H^EJy) - ^M^ilM 

k=l ^ 

I ^ --^^y^ - ^ -En-l{x) 



(1.9) 



^Bk{x)En-kiy) 
= i {Zl) ^ti|^B„_,(.)) (1.10) 



x-y J {x-yf 
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Corollary 1.3. Let n be any nonnegative integer. Then we have 

n+2 



(n + 2) 5^ E^{x)E^-u{x) -^Y.\^i W- l)-7^5n+2-/(x), 

fe=0 i=2 ^ ' ^ ' (1.11) 

^ ^^n-fc(a;) - T^-K^) = iiuEnix), (1.12) 

k=\ 1=2 ^ ^ 

n n / I 1 \ D 

5]Sfc(a;)£;,_fc(a;) - J] T ^ (2^ + Z - l)-^Er,-i{x) = {n+l)Er,{x). 

fc=0 /=2^^ + -^^ ' (1.13) 

Proo/. Letting y tend to x and noting that Ei{0) = 2{l-2^+^)Bi+i/{l + l), 
we then obtain the (1.11) and (1.12) from (1.8) and (1.9) respectively. 
Since 

En+i{y) - i?n+i(a:) = E'^+,{x){y - x) + ^^^^{y - xf + ■ ■ ■ , 



we have 



y^x \ X — y [x — yY 

En+i{y) - En+i{x) - (y - x)E'^+;^{x) 
■ ^™ ? \2 



2! 2 

Thus, (1.13) foUows from (1.10) by letting y tend to x. We are done. □ 
In the next section wc will prove Theorems 1.1 and 1.2 by a new ap- 
proach via differences and derivatives of polynomials. In Section 3 we will 
give a further extension of Corollary 1.1 with help of the gamma function. 

2. Proofs of Theorems 1.1 and 1.2 

Lemma 2.1. Let P{x), Q{x) G C[x] where C is the field of complex num- 
bers. 

(i) We have 

A{P{x)Q{x)) =P{x)A{Q{x)) + Q{x)A{P{x)) + A{P{x))A{Q{x)), 

=A*{P{x))A*{Q{x)) - P{x)A*{Q{x)) - Q{x)A*{P{x)) 

and 



A*{P{x)Q{x)) = A{P{x))A*{Q{x)) + P{x)A*{Q{x)) - Q{x)A{P{x)), 
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where the operators A and A* are given by A{f{x)) = f{x + 1) — f{x) and 



A*{f{x)) = f{x + l) + f{x). 

(ii) If A(P(x)) = A(Q(x)), then P'{x) = Q'{x). If A*{P{x)) = 
A*{Q{x)), then P{x) = Q{x). 



for every n = 1, 2, 3, . . . . Now that the polynomial g{x) = P{x) — Q{x) — 
P(0) + (5(0) has infinitely many zeroes, we must have g{x) — and hence 



P'ix) = Q'{x). 

Now assume that A*{P{x)) = A*{Q{x)). Then 

P{n) - Q{n) = -{P{n - 1) - Q(n - 1)) = . . . = (-l)^(P(O) - Q(0)) 



for every n = 1,2,.... Since the equations P{x) — Q{x) = P{0) — Q{0) 
and P{x) — Q{x) = — (P(0) — (5(0)) both have infinitely many roots, we 
must have P{x) — Q{x). 

The proof of Lemma 2.1 is now complete. □ 

Lemma 2.2. Let n be any positive integer. Then 



Proof, (i) Part (i) can be verified directly. 

(u) Suppose that A{P{x)) = A{Q{x)). Then, 



n— 1 n—l 



Pin) - P(0) = Yl ^(^(^)) = E ^(^(^)) = ^(^) - ^(0) 



k=0 i=0 




(2.1) 



and 



n 





'X 



.n—k 



.n—l 



(2.2) 



k=Q 



Proof Note that 
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where in the last step we apply a well-known identity of Chu (see, e.g. 
[GKP, (5.10)]). Similarly, we have 

fc=0 k=0 1=0 ^ ^ 

e-.(^)--t(:)=|:(::,>'(^)--' 



1=0 k=l ^ ^ 1=0 

So both (2.1) and (2.2) hold. □ 
Proof of Theorem 1.1. Observe that 



Therefore 



dxdyj^^\l-l) P 

d ^ fn-l\ Bi-i{y-x) 

= d^2:[l_l 1 Br.-i{x) 

1=1 ^ ^ 

E (^r^^)^5/-l(y-^)5n-l-/(^) 

E(^Ii^) ^5^-2(2/ -^)i?n-K^) 
1=2 ^ ^ 



1=2 

n-1 



1=2 
n-2 



1=0 



n-l\ Bi{y - x)Bn-2-i{x) 



l+lj 1+2 



d d ^(n-l\ Bi{x-y)Bn-i{y) 
dxdy^^\l-l) P 

d d ^ f n - 1\ Bi{x - y)Bn-i{y) 



dy dx ^ \l — 1 J P 

f n - 1\ Bi{x - y)Bn-2-i{y) 

1=0 



l + l J 1 + 2 
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We also have 



dx dy \ n{y — x) 
d ( nBn-i{y) Bn{y) - Bn{x] 



dx \ n{y — x) n{y — x)^ 

^ 2 B^{x)-B^{y) 



{y — x)"^ {y — x)'^ n {x — y)^ 

Let L{x,y) and R{x,y) denote the left hand side and right hand side of 
(1.4) respectively. In view of the above, that ^^L{x,y) = -^-^R{x,y) 
gives (1.5) with n replaced by n — 2. 

If we substitute a; + y for y in (1.4), then we get the following equivalent 
version of (1.4). 

Bk{x + y)Bn-k{x) 



k{n — k) 

= E^(;ri) my)Bn-i{x) + Bi{-y)Br,_r{x + y)) (2-3) 



k=i 

n 



n ny 

Now it suffices to prove (2.3) only. Let us view y as a fixed parameter. 
Denote by Pn{x) and Qn{x) the left hand side and the right hand side of 
(2.3) respectively. It is easy to check that 

p/ / ^ p / N Bn{x + y) + Bn{x) ^, / . ^ / ^ 

Pn+A^) - nPn[x) = = Qn+iix) - nQn{x). 

By Lemma 2.1(ii) we need only to show A{Pn+i{x)) — A{Qn+i{x)). 
In view of Lemma 2.1(i) and the fact that A{Bk{x)) = kx''~^, 

a(p„,.(x)) = |;a(^?^(^.|^) 



n-l 



fc=l fc=l fe=0 
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Applying Lemma 2.2 we then get 

[x + y) — X 

On the other hand, it is easy to see that A{Qn+i{x)) also equals the right 
hand side of the last equality. Therefore A(P„_|_i(a;)) = A{Qn+i{x)) as 
required. This concludes our proof. □ 

Proof of Theorem 1.2. Substituting a; + y for y in (1.8) we then get the 
following equivalent form of (1.8): 

i~ n Tb -t- y 

(2.4) 

'n + 1\ Ei{y)Bn+i-i{x) + Ei{-y)Bn+i-i{x + y) 



-2E 



I ) l + l 



1=0 

If we substitute x + y for x and x for j/ in (1.9), we then have the following 
equivalent version of (1.9): 

Bk{x + y) ^ , . , . ^n(a; + l/) -^n(a;) 

> , i En-k{X) - UnEn{X) 

^-^ K y 

(2.8) 

=11 (") (^^"-.W - ^4^^;„-,(x+.)). 

Note that (1.10) with n replaced by n — 1 follows from (1.9) by taking 
partial derivatives with respect to x. In view of the above, we only need 
to prove (2.4) and (2.5) with y fixed. 

(a) For ^ /c ^ n, by Lemma 2.1(i) we have 

/\{Ek{x + y)En-k{x)) 
=A*{Ek{x + y))A%En-k{x)) 

-Ekix + y)A*{En-k{x)) - En-k{x)A*{Eu{x + y)) 
=2{x + yf ■ 2a;"-^ - Ek{x + y) ■ 2a;"-^ - E,_fc(a;) • 2{x + yf. 

Thus 



Y,Ek{x + y)En-k{x)) -AY,{x + yf 

k=0 ' k=0 

n n 

= - 2 ^ Ek{x + y)a;-^ - 2 ^ En-k{x){x + y)\ 



fe=0 /c=0 
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With help of Lemma 2.2, we have 



n+1 _ rji.n+1 



From this we can easily check that A(P^(x)) = A{Qn{x)) where Pn{x) 
and Qn{x) denote the left hand side and the right hand side of (2.4) 
respectively. 

Now that A(P„+i(x)) = A((5„+i(x)), we have P^j^iix) = Q'^^^ix) by 
Lemma 2.1(ii). Clearly P^_^_i{x) = (n + 2)Pn{x) and Qn+i{x) — {n + 
'2)Qn{x)- So Pn{x) — Qn{x) and hence (2.4) holds. 

(b) Now we turn to prove (2.5). For 1 ^ /c ^ n, by Lemma 2.1(i) we 
have 

A*{Bkix + y)En-kix)) 
=A{Bk{x + y))A*{En-k{x)) 

+ Bk{x + y)A*{En-k{x)) - A{Bk{x + y))En-k{x) 
=k{x + yf-^2x''-^ + Bk{x + y)2x''-^ - k{x + yf-^En-k{x). 

Thus 

=2Y,{x + yf-'x-^ + 2 X: ^^Y^ ^"' " + y)'-'En-k{x) 

n-1 ^ X 
1=0 ^ ^ 

where we apply Lemma 2.2 in the last step. It follows that A*{L{x)) = 
A*(i?(x)) where L{x) and R{x^ are the left hand side and the right hand 
side of (2.5) respectively. Applying Lemma 2.1(ii) we find that L{x) = 
R{x). 

The proof of Theorem 1.2 is now complete. □ 
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3. Final remarks 

Quite recently Dunne and Schubert [DS] proposed the following gener- 
alization of Miki's and Matiyasevich's identities involving the well-known 
gamma function. 



i V 



^ B2kB2n-2k T{2k + p)T{2n-2k+p) 



r(2n + 2p) ^ 8k{n - k) T{2k)T{2n - 2k) 

-rl , 1^V^ B2kB2n-2k^i2k +p) B2n ^ST^ on , , ^\ 



2n-l 

(2A;)!(2n-2A;)!r(2/c + 2p+l) ' (2n)! E 
where n e {2, 3, . . . }, p ^ and 

Jo 



P{a,b) 



r(a + b) 



This result follows from our following theorem in the special case p = q 
and X = 0. 



Theorem 3.1. Let n > 1 be an integer, and let p ^ and q ^ 0. Then 

Tin + p + q) 

^fn-l\B, r(l+p)T(q+l)+T(l + q)T(p+l) (3.1^ 

= E(^,_ijT^n-K-) r(/ + p + , + i) 



where 



n—l f 

Hn{p.q) = Y^(3{k + q,p+l) = \ 
k=l ^ 



P{p,q + 1) - (3{p,n + q) ifp^O, 



In the case p — q — Ooip — q — 1 Theorem 3.1 yields (1.6) or (1.7). 

We omit our proof of Theorem 3.1 which is very similar to the proof of 
Theorem 1.1. The key point is the following extension of Chu's identity 
which can be easily proved by induction on n. 

Lemma 3.1. Let n ^ I > be integers, and let p ^ and q > 0. Then 
E Z J) Pik + P,n-k + q) = (3{l +p,q). (3.2) 
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